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ABSTRACT 

Context. Broad absorption line quasars (BAL QSOs) belong to a class of objects not well-understood as yet. Their UV spectra show 
BALs in the blue wings of the UV resonance lines, owing to ionized gas with outflow velocities up to 0.2 c. They can have radio 
emission that is diificult to characterize and that needs to be studied at various wavelengths and resolutions. 
Aims. We aim to study the pc-scale properties of their synchrotron emission and, in particular, to determine their core properties. 
Methods. We performed observations in the Very Long Baseline Interferometiy (VLBI) technique, using both the European VLBI 
Network (EVN) at 5 GHz, and the Very Long Baseline Array (VLBA) at 5 and 8.4 GHz to map the pc-scale structure of the brightest 
radio-loud objects of our sample, allowing a proper morphological interpretation. 

Results. A variety of morphologies have been found: 9 BAL QSOs on a total of 1 1 observed sources have a resolved structure. 
Core-jet, double, and symmetric objects are present, suggesting different orientations. In some cases the sources can be young GPS 
or CSS. The projected linear size of the sources, also considering observations from our previous work for the same objects, can vary 
from tens of pc to hundreds of kpc. In some cases, a diifuse emission can be supposed from the missing flux-density with respect to 
previous lower resolution observations. Finally, the magnetic field strength does not significantly dilfer from the values found in the 
literature for radio sources with similar sizes. 

Conclusions. These results are not easily interpreted with the youth scenario for BAL QSOs, in which they are generally compact 
objects still expelling a dust cocoon. The variety of orientations, morphologies, and extensions found are presumably related to 
different possible angles for the BAL producing outflows, with respect to the jet axis. Moreover, the phenomenon could be present in 
various phases of the QSO evolution. 

Key words. Quasars: absorption lines - Galaxies: active - Galaxies: evolution - Radio continuum: galaxies 



■1. Introduction 



^ The nature and origin of broad absorption line quasars (BAL 
• ^H QSOs) is still an open issue in the framework of active galac- 
^^ tic nuclei morphology and evolution. The key characteristic 
H of this class of objec ts (~20% of the entire QSO population, 
i . . jHewett & Foltzll2003h resides in their UV spectra, where broad 
absorption lines (BALs) are present in the blue wings of the U V 
resonance lines, owing to ionized gas with outflow velocities up 
to 0.2 c. Some evidence of BALs in young or lumin ous IR galax- 
ies has led to propos ing an evolutionary scenario (iBriggs et alJ 
Il984t ISandersll2002h . in which this phenomenon is due to the 
young age of these objects: the central AGN would still be ex- 
pelling the enveloping dust cocoon. The fraction of the QSOs 
showing BALs can then be interpreted as the duration of this 
phase in relation to the total active life. This view was supported 
by the presence of prototypical characteristics of young radio 
sources, like GigaHertz Peaked Sourc es (GPS) or compact steep 
sources (CSS), among BAL QSOs (iMontenegro-Montes et al.l 
l2008al:lLiuetal.ll2()081) . 



An orientation model proposed by |Elvisl(l2000l) predicts that 
BAL outflows are present in all QSOs, but only when they inter- 
cept the line of sight would the BALs be detected. In this case the 
percentage of BAL QSOs would constrain the solid angle of t he 
outflow. The bipolar wind model pro posed by lPunslvl (Il999al[bh 
was used by iGhosh & Punslv (l2007h to explain the polar BAL 
QSOs found by lZhou et alj (l2006h :in this model the BAL out- 



flow would be aligned with the polar axis, while the relativistic 
jet would be nested inside of it. Recent works, based on large 
samples, do not confirm a young age for all BAL QSOs, but 
also a variety of orient ations have been found (iBruniet al.ll2012l : 
iDiPompeo et alJl201 ih . Also, on the arcsec-scale the majority 
of these sources remain unresolved in the radio ban d, with only 
about 10% or less showing an extended morphology (IBruni et al.l 
120 121: IDiPompeo et alJl2()llh . 

Radio-loud BA L QSOs are a sniaU fra c tion of the BAL QSO 
population (-10%. iHewett & Foltzl 120031: IShankar et al.ll2008l) . 
but the radio emission can be a n additional tool for un derstand- 
ing the orientation and the age (lO'Dea & Baum|[l997l) of these 
sources. We focus on this class of objects to collect information 
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that can help in understanding which of the two models is the 
most probable. Thus, we embarked on a VLBI project to test the 
inner structure of the QSO. 

The VLBI technique is an important tool for this kind of 
study, allowing one to reach the necessary resolution and sen- 
sitivity to study distant objects like BAL QSOs. This approach 
offers different indicators of the morphology of the radio source: 
(1) The detection of radio jets and their asymmetry can provide 
constraints on the radio source orientation; (2) The spectral in- 
dex analysis allows us to discriminate between core-jet or dou- 
ble structure; (3) Variability with respect to the total flux-density 
already measured with the VLA can give an estimate of the ac- 
tivity and orientation of the QSO jet axis with respect to the line 
of sight. (4) In case the pc-scale structure displays mini-lobes, 
which are typically found in young radio sources, the projected 
linear size of the radio source ca n be used to estimate the kine- 
matic age (iDallacasa et alj|2003h . 

In previou s works on BAL QSOs, using the VLBI imag- 
ing technique d Jiang & Wan a 2003 ;_ Kunert-Bairaszewska et al.l 



2007*. '20 101; iLiu et al.1 12 008'; 'Montenegro- Montes et al.ll2008bl; 
Gawrohs ki & Kunert-Bai raszewska 2011), various morpholo- 
gies and sizes have been found, and in most cases sources were 
unresolved even at pc-scale resolution. In this paper we present 
the first results of an extensive observational campaign to deter- 
mine the pc-scale structure and morphology of a complete sam- 
ple of radio-loud BAL Q SOs, whose characteristics in the radio- 
band were presented in iBruni et all (l20I2b . Beyond the study 
of the pc-scale radio properties of the BAL QSO population, we 
aim at distinguishing which of the presented models could be the 
most suitable for explaining the BAL phenomenon. 

This paper is organized as follows. In Sec. 2 we present the 
sources and the observations, and in Sec. 3 we show the pc-scale 
radio maps and discuss the morphology. In Sec. 4 we discuss our 
results in the light of earlier works in the literature. 

The cosmology adopted throughout the paper assumes a flat 
universe and the following parameters: Hq-IQ km s ' Mpc"', 
f2A=0.7, Qm=0.3. The adopted convention for the spectral index 
definition is Sv oc v". 



2. Radio observations 

This work is based on a well-defined s ample of radio-loud BAL 
QSOs presented in iBruni et aP (l2012l) . This sample has been 
selected from the fourth rele ase of the Sloan Digiti zed Sky Sur- 
vey (SDSS) QSO catalogue JSchneider et al.ll2007l) . with a flux- 
density threshold of 30 mJy at 1.4 GHz and a redshift range 
of 1.7<z<4.7: these criteria resulted in a sa mple of 25 objects . 
The radio spectra of the sources, presented in lBruni et al.l(l2012l) . 
show that BAL QSOs are weaker at radio-frequencies above 1 .4 
GHz, with a few exceptions, requiring sensitive observations to 
study their pc-scale structure. 

In this paper, we present VLBI observations of 11 of the 
brighter sources [S 1.4 >~ 50 mJy) included in our sample. 
These 1 1 are a random subset of these brighter sources. The 
high-resolution study of sources 1159+01 (5 1.4 = 266 mJy) 
and 1624+37 (5 1. 4 - 56 mJy), not included here, has already 
been presented by iMontenegro-Montes et al.l (|2008b): we thus 
obtained parsec-scale radio maps for a total of 13 out of 25 BAL 
QSOs in the whole sample. The study of the polarized flux den- 
sity was not possible in these works, since much more observing 
time would have been needed to detect the low fraction of polar- 
ized flux density of these objects. More information on overall 
radio spectra and the kpc-scale polarization properties of the ob- 
jects studied here, as well as the remaining sources in the sample. 



are available in iBruni et al.l (l2012l) . In Table [T] we report some 
basic information on the 1 1 BAL QSOs presented here, while a 
summary of the observations is reported in Table |2] 

All the data reduction and map analysis was performed using 
the NRAO AIP^ software, following the usual procedures once 
editing of raw visibilities and of antenna temperature (TSYS) 
measurements were carried out. Based on the variation in the 
amplitude solutions once system temperatures and antenna gain 
have been considered, the accuracy of the flux-density scale can 
be estimated to be within +5%. The final angular resolution of 
the maps is in the range 2-5 mas. 

2.1. VLBAdata 

Seven objects were observed with the Very Long Baseline Ar- 
ray (VLBA) at 5 and 8.4 GHz (C and X-band, respectively) with 
a total bandwidth of 64 MHz at each frequency, in two observ- 
ing runs in 2010 and 201 1. In the first run an average on-source 
time of about one hour at 5 GHz and two hours at 8.4 GHz was 
allowed. The correlation was performed with the VLBA cor- 
relator at the National Radio Astronomy Observatory (NRAO) 
in Socorro (US). In the second run, when the fainter sources 
were observed, a total of three hours at 5 GHz and four hours 
at 8.4 GHz were spent on each target. The data were then pro- 
cess ed with the new D istributed FX (DiFX) software correla- 
tor (iDelleretal .11201 ll) . Given that, in general, the target BAL 
QSOs were weak; the VLBA observations were carried out in 
phase referencing mode for most of sources; the transfer of de- 
lay, phase, and phase rate solutions from a reference source close 
to target allowed absolute positions to be obtained, which signif- 
icantly improved the earlier infor mation taken from the FIRST 
catalogue (Becker et al. 1995; Whi te et al.l 1997h . Therefore, the 
phase centre in the maps presented in Figure prefers to the shift 
applied and not to the peak in the map plane whose position is 
instead reported in Table [1] Since the correlation was done us- 
ing the position from the FIRST survey, we had to allow a rather 
wide field for the initial imaging of each target source. In fact 
we started with about 1" wide maps, and once the target source 
was detected, a smaller field was selected after applying appro- 
priate shifts in RA and DEC in order to have the target source at 
about the centre of the map. The same shift was applied at both 
frequencies. Then a few iterations of phase self-calibration were 
performed. The final maps are shown in Fig. |3]and are discussed 
in the next section. 



2.2. EVNdata 

Four additional sources were observed with the European VLBI 
Network (EVN) at 5 GHz (C-band), in 2009, using the anten- 
nas in Jodrell Bank, Westerbork, Effelsberg, Onsala, Medicina, 
Torun, Shanghai, Urumqi, Noto, Yebes, and the two MERLIN 
telescopes in Damhal and Knockin. A total bandwidth of 32 
MHz and an average on-source time of about two hours were al- 
lowed. The data were correlated at the Joint Institute for VLBI 
in Europe (JIVE) correlator in Dwingeloo (The Netherlands). 

For three of the sources we carried out phase-referencing. 
The maps thus adopt the same convention as for the VLBA ones. 
Standard procedures, including fringe-fitting, were used for the 
data reduction of the target sources. Delay, phase, and rate so- 
lutions were generally found above the standard signal-to-noise 
ratio threshold, for all the target sources. 



http://www.aips.nrao.edu 
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Table 1: Sample of the 1 1 radio-loud BAL QSOs studied in this paper. 



Name 


SDSS name 


RA 


DEC 


z 


5 1.4 BAL 






(J2000) 


(J2000) 




(mJy/beam) Type 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) (7) 



0044-hOO 


J004444.06+001303.5 


00 44 44.0657 


+00 13 03.560 


2.28 


53.1 


- 


0756-H37 


J075628.25+371455.6 


07 56 28.24 


+37 14 55.6 


2.51 


239.4 


- 


0816-H48 


J081618.99+482328.4 


08 16 18.9960 


+48 23 28.469 


3.57 


68.3 


- 


0849+27 


J084914.27+275729.7 


08 49 14.2716 


+27 57 29.709 


1.73 


52.8 


Hi 


1014+05 


J101440.35+053712.6 


10 14 40.3561 


+05 37 12.682 


2.01 


55.0 


Hi 


1102+11 


Jl 10206.66+112104.9 


1102 06.6617 


+ 112104.900 


2.35 


82.3 


- 


1237+47 


J123717.44+470807.0 


12 37 17.44 


+47 08 07.0 


2.27 


78.5 


FeLo 


1304+13 


J130448.06+130416.6 


13 04 48.0594 


+ 13 04 16.590 


2.57 


49.6 


- 


1327+03 


J132703.21+031311.2 


13 27 03.21 


+03 13 11.2 


2.83 


60.7 


- 


1406+34 


J140653. 84+343337.3 


14 06 53.84 


+34 33 37.3 


2.56 


164.4 


- 


1603+30 


J160354. 15+300208.6 


16 03 54.1534 


+30 02 08.705 


2.03 


53.7 


Hi 



Notes: Column 3,4, and 5 are respectively RA, DEC and redshift as measured from the SDSS. For phase-referenced sources we 
give the coordinates of the brightest compon ent peak, derived from our VLBI ob servations. Column 6 gives the peak flux 
d ensity at 1.4 GHz fr om the FIRST catalogue dBecker et al.lll995l:rWhite et alJll997h . The last column is the BAL type as given 
in lBruni et al.l ( |2012|) : a hyphen indicates that it could not be determined. 



Table 2: Summary of the observations and instrumental setup, 
power beam-width). 



In the penultimate column the typical beam sizes are given (half- 



Run 


Date 




Telescope 


Frequencies 
(GHz) 


Bandwidth 
(MHz) 


Ghpbw 
(mas) 


Time on source 
(hour) 


1 

2 
3 


25-26 Oct 2009 

19, 22, 26, 28 Feb 2010 

4, 7, 23, 26 Apr -1,2 May 


2011 


EVN 
VLBA 
VLBA 


5 
5,8.4 
5,8.4 


32 
64 
64 


5 
4,2 
4,2 


2 

1,2 
3,4 


ilts 








0044+00 









Once final maps had been obtained, we carried out a 2-D Gaus- 
sian fitting solving for component position, flux-density, size, 
and position angle. The resulting information is reported in Ta- 
ble [3] The distance between the centroids of each individual 
component as found by the fits were used to calculate the total 
projected angular and linear sizes of the source. 

In the following, we briefly describe each source, observed 
with the VLBA (Section 3.1) or the EVN (Section 3.2). For the 
often mentioned overall radio-spectrum of the sources we re- 
fer to lBruni et al.l (|2012|) . Figure [T]provides a comparison with 
the VLA flux-densities presented in our previous work. In this 
section we also discuss any hint of variability for the individual 
sources. Finally, in Table[T]we provide the position of the bright- 
est component peak for the seven sources with absolute position 
information from phase-referencing. 



3.1. VLBA maps 

In this section, the pc-scale structure of the BAL QSOs observed 
with the VLBA is briefly discussed. Most of the sources have 
been found to be unresolved in the earlier arcsecond-scale ob- 
servations. The typical resolution of the VLBA maps is ~ 4 mas 
at 5 GHz and ~ 2 mas at 8.4 GHz. The spectral indices af* for 
the individual components were calculated using the integrated 
flux density from the Gaussian fit, and do not refer any to local 
values. The availability of maps at two frequencies often helps 
in correctly interpreting the source morphology in terms of cores 
or steep spectrum components. 



Both maps at 5 GHz and 8.4 GHz show a resolved structure. In 
particular at 5 GHz, three misaligned components are present: 
two of them are clearly visible, the third is embedded in some 
diffuse structure. While a Gaussian fit was possible for com- 
ponents A and B, for component C it did not converge since it 
barely stands against some diffiise emission, and only an upper 
limit for its flux-density could be obtained. In case the map at 
8.4 GHz was not available, component C would not have been 
considered. In this 8.4 GHz map, much deeper than that at 5 
GHz, only components A and C are visible, suggesting a steep 
spectral index for component B whose surface brightness drops 
below the detection level. Component A has a rather flat spectral 
index (-0.52+0.27) and could be most probably classified as the 
core region, in which part of the jet is also visible, as suggested 
by the elongation in p. a. ~ 176°. 

Although it is not straightforward, the morphology of this 
source can be interpreted as a core-jet. 

0756+37 

In both VLBA bands we resolve the radio source into two com- 
ponents, separated by ~1.8 mas (~15 pc at the redshift of the 
object) at the highest resolution, in p. a. ~ -60°. Both compo- 
nents have similar sizes and spectral indexes (-0.38+0.27 for A 
and -0.37+0.27 for B). Since the overall spectrum peaks at 2.5 
GHz (8.7 GHz rest frame), and it is steep in the optically thin 
(power-law, not self-absorbed) part of the emission, we can in- 
fer that we have detected two hot spots in the mini-lobes of a 
typically young radio source. 
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Fig. 1 : Flux densities of VLBI vs VLA at 5 (left) and 8 GHz (right), for all of the sources with both measurements. The dashed line 
indicates a slope equal to 1. The outUers {(T,,ar > 3) are sources 0756+37, 1406+34 and 1014+05 at 5 GHz, and 1237+47 at 8.4 
GHz. 



The VLBA can account for a fraction of the total flux den- 
sity measured in earlier VLA observations (~65% at 5 GHz and 
~85% at 8 .4 GHz): c o nsider ing the variability significance as 
defined by IZhou et al] ( |2006|) , the discrepancy is significant at 
5 GHz, with a (Xvar^T-O, while at 8 GHz we obtain cr,,^,,. < 3. 
The time interval between VLA and VLBA observations is about 
seven months. Thus, either the source is variable - but this would 
be inconsistent with the hot-spot scenario - or, more likely, there 
is some emission on angular scales not sampled by the VLBA 
(tens of mas or more). 

This source was o bserve d also by 

iGawronski & Kunert-Bai raszewskal (1201 ih using the EVN 
at 1.6 GHz, and reporting a probable core-jet structure, but no 
maps were provided in that work. 

0849+27 

This source has a lready been resolve d in the FIRST survey (VLA 
at 1 .4 GHz, see iBruni et al] 1201 2h showing three components 
with a projected linear size of 382 kpc, and a position angle for 
the structure of ~ 45°. These VLBA observations were designed 
to study the central component, which was unresolved on the 
arcsecond scale. The maps show two components at 5 GHz (A 
and B) separated by 3.5 mas (30 pc), while only B is visible at 
8.4 GHz, with a spectral index of -0.12+0.20 (flat) that can be 
considered as the source core. This component is resolved at 
both frequencies in p. a. ~ 30°. 

1102+11 

The map at 5 GHz shows a structure resolved into two compo- 
nents, which are separated by -4.1 mas (~34 pc), in p. a. -53°. 
At 8.4 GHz only the strongest component (A) has been detected, 
and it turned out to be resolved along the same direction. Its 
spectral index is -0.97+0.29. Such a steep spectrum makes it 
unlikely that component A is the source core. It is also too steep 
to be an active hot spot. Component B has probably an even 
steeper spectral index, since the flux density is below the 3-cr 
significance limit at 8.4 GHz. All this makes the morphological 



classification of the pc-scale structure rather unclear 
The source presents a convex radio spectrum, although more ob- 
servations at low frequencies are needed to clearly determine the 
turnover frequency, if present. 

1237+47 

This source is resolved into three components at 5 GHz, with 
the strongest one in the centre (A). At 8.4 GHz, the eastern 
component C is resolved in two subcomponents (CI and C2), 
both aligned with A and B. Component A has an inverted 
spectral index (0.41+0.29) and thus could be interpreted as the 
core, while B and CI present elongated emission toward A, with 
a similar flux density at 8.4 GHz. It could be the inner part of 
the jets, with an additional hot spot (C2) visible in the direction 
of CI. The projected linear size from X band is 4.8 mas (40 pc). 
From our previous radio observations of this source we know 
that it presents a flat spectrum, which can possibly be approxi- 
mated with a convex shape with a peak at about 3.5 GHz. Thus, 
this object can be a young radio source in which the core region 
is still the brightest in the source. Moreover, at 8.4 GHz its flux 
density exceeds (~37%) the total flux density in our earlier VLA 
observation. This could indicate that there is still a very active 
phase in which new components are likely to be formed on 
short timescales, creating the variability in the flux density. The 
time interval between VLA and VLBA observations is about 
seven months, with a cr„,^=3.4, just above the 3-cr significance 
threshold. 



1304+13 

It was possible to resolve this source both at 5 GHz and at 8.4 

GHz: in the maps, the two brightest components are on the edges 

of the radio source (A and B), while the central component C is 

distinguishable only at 8.4 GHz. The fit did not converge for this 

last component, so we could only estimate the flux density upper 

limit. 

The total projected linear size of the source at 8.4 GHz is 5.3 
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Fig. 2: Maps of the 6 BAL QSOs resolved with the VLBA. The synthesized beam size is shown in the lower left comer of the map. 
Levels are 3-cr multiples, according to the legend. Orientation is N-up, E-left. 
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Fig. 2: Continued. 
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mas (44 pc). The overall spectrum is steep, and does not present 
a peak in the GHz range, which is the only one that could be 
investigated at the moment. 

1406+34 

This source is unresolved in both bands with the VLBA. We 
can set an upper limit of 1.8 mas (14 pc) to its size. The map 
at 8.4 GHz shows a noise pattern that is not related to any ex- 
tended structure. The pc-scale emission revealed by the present 
observations has an inverted spectrum (q'=0. 36+0.26), and we 
are therefore seeing the optically thick, self-absorbed part of the 
synchrotron emission. 

The overall radio spectrum peaks at about 5.4 GHz. The 
VLBA flux-density measurements at 5 GHz accounts only for 
~67% that of the VLA (cr,,Q^=4.8), owing either to some intrin- 
sic variability or to an additional component fully resolved on 
the pc-scale. Such a component, with a steep spectrum, may be 
respons ible for the 327 MHz emission detected in the WENSS 
survey dde Bruvn et al.l 120001) . which could be the remnant of 
an earlier episode of radio activity and which also modifies the 
spectral slope at low frequencies. The source is also unresolved 
in VLA maps. 

3.2. Results from EVN observations 

The EVN observed four BAL QSOs at 5 GHz (C-band). This 
array is made of various radio telescopes with substantially dif- 
ferent performances. In particular, the presence of a large dish 
like the Eff'elsberg-lOOm allows better baseline sensitivity than 
on a typical VLBA baseline. However, due to the geographi- 
cal distribution of the telescopes, the uv-plane is not uniformly 
sampled, and, on average, the structural information that can be 
derived here is poorer than the maps presented in the previous 
section. The typical resolution of these maps is ~5 mas. 

0816+48 

This source shows an elongation toward SE. From the major axis 
of the Gaussian fit (6.8 mas) we can infer an upper limit of pro- 
jected linear size of 51 pc. A clear classification of the spectrum 
was not possible in our previous work, but the VLA map at 1 .4 
GHz already showed an elongation toward SW, not correspond- 
ing with the one present here, and with a major and minor axes 
of 29 and 14 arcsec, corresponding to 217 and 105 kpc, respec- 
tively. This could be another case of intermittent activity with 
a change in the jet orientation between the two phases. Further 
observations at the intermediate VLA resolution could confirm 
this hypothesis. 

1014+05 

This source presents a complex morphology whose proper clas- 
sification is rather difficult, because the brightest component (A) 
can be interpreted as the core region, and there are two more 
components (B and C) toward north. The overall structure is 
well aligned in p. a. ~0°. 

The distance between the centroids of A and C components 
is ~26.6 mas, corresponding to a projected linear size of 229 
pc. In our earlier work we found a power-law spectrum for this 
source, but more data below 1 .4 GHz are needed to exclude any 
turnover in the MHz range. 
Our EVN map can account for ~59% of the total flux density 
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measured in the VLA map (cr,,a,.=8.2), since this discrepancy 
may arise either from some flux-density variability or from an 
additional structure not sampled by the EVN spacings. 

1327+03 

This is an unresolved source. The upper limit for the linear size, 
estimated from the Gaussian fit, is 46 pc (5.9 mas) and 11 pc 
(1.4 mas), respectively, on the major and minor axes, but the 
beam is considerably elliptical in this map so this can only be an 
indicative value. A power-law spectrum down to 1 .4 GHz does 
not show any turnover for this source. 

1603+30 

In this case a rather symmetric structure can be seen on the 
parsec scale. The south-western component can be fitted with 
two (A and B) Gaussians. The separation between A and C (at 
the north-eastern edge) components is 8.6 mas, with a projected 
linear size of 74 pc. 

This source presents a convex radio spectrum, with a peak 
frequency at 1 .4 GHz, and thus could be classified as a young 
GPS source, in which the emission is dominated by the hot spots. 

On the arcsecond scale, a VLA map at 22 GHz shows a re- 
solved structure in which some additional emission is present at 
about 2 arcsec (17 kpc) to the south. Such a linear size is more 
typical of CSS sources (1-20 kpc), although the overall spectrum 
is typical of GPS sources. The structure seen with the EVN is 
not aligned with the one seen with the VLA; probably a merger 
or a jet precession could justify the different morphologies of the 
same object on different scales. This is another case of a possi- 
bly young source, like 1237-1-47, with a flux density at VLBI 
resolution exceeding the earlier VLA observations at the same 
frequency (~25%), but in this case the variability significance 
(Tvar is Only 2.5, i.e. below the threshold. This sour ce has al- 
ready been observed with EVN+MERLIN at 1.6 GHz ( iLiu et alj 
|20085), but owing to the lower resolution (13.0 x 4.0 mas) it was 
unresolved. 
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Table 3: Results for the 1 1 radio-loud BAL QSOs studied in this paper. 



Name comp. 



(1) 



(2) 



gVLBI 

(mJy) 
(3) 



oVLBl 
■^8.4 

(mJy) 
(4) 



(5) 



cVLA 
''4.8 

(mJy) 
(6) 



cVLA 
'^84 

(mJy) 

(7) 



Maj. axis Min. axis Maj. axis Min. axis P. A. LS LS 

(mas) (mas) (pc) (pc) (degrees) (mas) (pc) 

(8) (9) (10) (11) (12) (13) (14) 



S«, 



B'eq 



Morphology ((„„■„ ^«, 

(lO-^'erg/cm^) (mG) (mG) 
(15) (16) (17) (18) 



> 
o' 
B 






0044+00 



0756+37 



0816+48 
0849+27 



1014+05 



1102+11 



1237+47 



1304+13 



1327+03 
1406+34 
1603+30 



A 
B 
C 

A 
B 



A 
B 

A 
B 
C 

A 
B 

A 
B 
C 

CI 
C2 

A 
B 
C 



A 
B 
C 



33.0+2.6 
25.0+2.5 

6.5+0.8 
<2.4 

147+11 
88.2+8.8 
59.2+5.9 
26.9+2.7 
14.8+1.3 

2.2+0.3 
12.6+1.3 
20.4+1.6 
10.3+1.1 

7.4+1.0 

2.7+0.6 
38.8+3.4 
28.1+3.0 
10.7+1.7 
54.4+4.5 
37.8+4.0 
11.0+1.5 
11.4+1.7 



20.1 + 1.6 

7.2+0.7 

12.9+1.4 

80.6+8.0 
211+21 

43.3+2.6 

18.1 + 1.8 
8.8+1.0 

16.4+1.6 



20.8+1.9 
19.0+1.9 

1.8+0.2 
121+9 

72.2+7.2 
48.8+4.8 

11.8+0.2 

<0.6 
11.8+0.2 



16.9+1.8 
16.9+1.8 

<5.4 
84.6+6.6 
46.9+5.0 
16.1+2.6 

14.3+2.9 
7.3+1.7 

15.1+0.5 
3.2+0.3 

11.2+0.3 
<0.9 



34.2+1.8 25.5+0.6 



-0.52+0.27 

>-0.5 

-0.38+0.27 
-0.37+0.27 



255+25 0.36+0.26 



226+2 142+2 



31.4+0.5 19.5+0.5 



<-2.5 
-0.12+0.20 



34.5+0.6 25.1+0.5 



39.8+0.8 19.9+0.7 



-0.97+0.29 
<-1.3 

0.41+0.29 
0.72+0.40 



62.3+1.0 61.6+1.6 



24.6+0.5 17.1+0.5 



-1.54+0.26 
-0.27+0.21 



79.5+1.7 56.5+0.7 
313+3 276+3 
34.5+2.4 26.9+0.6 



1.3 
4.5 
2.9 

0.3 
0.6 

3.3 

<2.9 

2.7 

<7.3 
6.7 
2.9 

<2.8 
2.4 

1.2 
1.8 

2.6 

2.2 

0.8 
1.1 

<5.9 
<1.8 

2.4 
2.4 

2.2 



0.7 
1.5 
0.7 

0.1 
0.5 
0.1 

<1.1 

0.4 

<5.3 
2.6 
0.9 

<0.9 

1.2 

0.8 
0.6 

1.6 

0.3 

0.5 
0.3 

<1.4 
<0.7 

1.3 
1.4 
1.5 



11 

37 
24 

2 

5 

24 

<25 

23 

<61 
56 

24 

<23 
20 

10 
15 

21 
18 

6 
9 

<46 
<14 

20 
20 
18 



6 

13 

6 

1 
3 
1 

<9 

3 

<45 

22 



<7 
10 

7 
5 

13 

2 

4 

2 

<11 
<6 

11 
12 
13 



176 

8 

171 

6 

8 
108 

26 

48 

155 
165 

122 

167 
10 

5 
178 

11 

179 

3 
174 

9 

5 

18 
18 
15 



5.6 47 



15 



6.8 

3.5 



4.1 



<5.9 

<1.8 



51 
30 



26.6 229 



34 



4.8 40 



5.3 44 



<46 
<14 

74 



CJ 



D 



CJ 
CJ 



M 



D 



M 



CJ 



u 
u 

M 



1.1 


3.4 


2.4 




0.4 


2.0 


1.5 




1.1 


3.4 


2.4 




7.7 


9.1 


5.7 





2.6 


5.3 


3.5 


ffl 


10.5 


10.6 


6.5 


3^ 


>0.4 


>2.0 


>1.5 




1.4 


3.9 


2.7 


H 


>0.1 


>1.0 


>0.8 


13 


0.2 


1.4 


1.1 




0.5 


2.3 


1.7 




£5 


>1.0 


>3.3 


>2.3 




0.5 


2.4 


1.7 





2.0 
4.0 

>0.6 



0.4 
0.3 
0.3 



4.6 
6.6 



3.2 
4.3 



>2.6 >1.9 



2.0 
1.9 
1.8 



1.5 
1.5 
1.4 



Notes: In column 2 the component is specified, colum ns 3 and 4 are flux densities measured with the VLBA or EVN; column 5 is the spectral index for the components; columns 
6 and 7 are the VLA flux densities as presented in lBruni et al.l(l2012l) : columns 8-11 are the deconvolved major and minor axes of the Gaussian fit, at the highest resolution 
map available, for each component. Column 12, 13, and 14 are the position angle of the single components (from the same map used to calculate the component size) and 
the total projected linear size of the source. Column 15 is the morphology classification (CJ; core-jet; D; double; M; multiple; U; unresolved). Columns 16, 17, and 18 are 
the derived quantities from the magnetic field estimation (see formulae [T]|2] and|3]l. 
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3.3. Magnetic field estimation 



We applied the iGovoni & Ferettil ( |2004|) formulae to calculate 
the magnetic field intensity or strength B. Supposing the syn- 
chrotron mechanism to be at the origin of the radio emission, 
this estimate assumes that the magnetic-field energy density Ub 
and the relativistic particles energy density Upari, contained in 
the emitting region, contribute with approximately equal values 
to the total energy density t/,o,. This equipartition assumption 
allows minimizing [/,„,, providing a conservative estimate of B. 
The formulae for minimum energy density, m,„„,, and the corre- 
sponding equipartition value for the magnetic field strength, Beq, 
are the following: 



UminV^] = ^(a,Vi,V2)(l +^r'(yO[MHz]r"(l + z) 



,4/7^,, „ ,^4ff/7^j ^ ,s^(i2+4a)/7 



x(/o[^i_]r'(4^pc]r 



B. 



_ 2An 

eq — v~^^ ^n 



1/2 



(2) 



where z is the source redshift, /o the source brightness at the 
frequency vq, d the source depth, k the proton-to-electron energy 
ratio, vo the observed frequency, and a the spectral index of the 
energy distribution of the electrons. The constant ^(a, vi, V2), 
depending on a and on the minimum and maximum energy of 
the charged particles (corresponding to v\ and V2), is tabulated 
in Table 1 of loovoni & Ferettil (120041) . We chose a value of 2.5 x 
10"'^ corresponding to a = 0.5, and minimum and maximum 
emission frequencies of 10 MHz and 100 GHz, respectively. In 
this calculation we assumed A;=l, a magnetic field filling-factor 
0=1, and an ellipsoidal volume for the component, with depth d 
equal to the minor axis. 

Assuming a low-energy cut-off in the particle energy distri- 
bution, the derived magnetic field B' is related to Bgq by the 
following formula (iBrunetti et al.ll 19971) : 



"eq 'mm "eq ■ 



(3) 

We assumed a minimum Lorentz factor 7,„,„ = 100 in this calcu- 
lation. 

In Table [3] we present the results of the magnetic field esti- 
mate for each component (or for the all the sources in case it is 
unresolved). For the unresolved sources or components, having 
deconvolved dimensions equal to zero, we used the major and 
minor axes of the beam and the peak flux density. Components 
without Gaussian fit were omitted. For two sources (1237+47 
and 1406-1-34), the estimation was not possible, since compo- 
nents show an inverted spectrum, which means we are looking at 
the self-absorbed part of the spectrum: in this case, the measured 
flux density is biased by the absorption, and does not reflect the 
real emission of the component. 

We found values of m„„„ ranging from 10"^ to 10"^ erg/cm-', 
corresponding to B' between 0.8 and 6.5 mG. In Sec. 4 we 
compare these results with values from literature. It is worth 
noting that some a priori assumptions have been made in this 
estimation. The value of k, the filling factor of the magnetic field 
(p, and the depth of the source along the line of sight d can have 
slightly diff'erent values from the ones used here, changing the 
results up to an order of magnitude. 

4. Discussion 

Radio-loud BAL QSOs were considered extremely rare sources 
before the advent of large surveys like the SDSS. After the re- 
lease of SDSS DR4, it has been possible to select larger samples 
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Fig. 4: Integrated flux density from the FIRST versus red- 
shift for the 13 BAL QSOs with VLBI observations com- 



prised in the Bruni et al. (|2012 ) sample (from this paper or from 
iMontenegro-Montes et aljr2008bl - crosses) and from the litera- 
ture (circles). For source 0849-1-27, resolved at the FIRST res- 
olution, we considered the flux density given by the sum of the 
single components. 



of BAL QSOs with a considerable emission also in the radio do- 
main. The sample we studied here has been defined to have the 
largest possible number of sources with radio flux density high 
enough to allow a detailed study of the pc-scale structure. 

In the literature, 32 BAL QSOs have to date 
been observe d with the VLBI imaging t ec hnique 
dJiang & Wand 120031: |Kunert-Bairaszewska et al.' '2007i 



2010t iLiuetaLll 



iKui 
jMc 



2008 | : | Montenegro-M ontes et al. 2008b; 

Gawronski & Kunert-Bairaszewskall201 lb . With this work we 
add a further nine previously unobserved sources, improving 
the statistics concerning the pc-scale radio morphology. Only 
8 BAL QSOs out of 32 (~25%) were found to be unresolved 
in earlier works, while the remaining objects show a variety of 
morphologies, including core-jet, quite symmetric, or complex 
structures. In some cases, a re-orientation or jet-precession 
scena ri o has been proposed (iKunert-Bairaszewska et al.l 120071 
I2OIOI : iGawroriski & Kunert-B airaszewskal | 2pllh to explain 



the pecuUar morphology. iDoi et al. ( 20091) have observed 



another 22 BAL QSOs with the Optically ConnecTed Array for 
VLBI Exploration (OCTAVE), subarray of the Japanese VLBI 
Network, to detect pc-scale radio emission due to non-thermal 
jets. Their work does not provide maps, hence no morphological 
information comparable to the previous, but also in this case 
the results did not point to a single pole-on geometry for BAL 
QSOs. 

Concerning our sample, four sources (16%) out of the 25 in 
the whole sample (0816+48, 0849+27, 1103+11, and 1603+30) 
were already resolved at the inte rmediate angular resolution of 
our previous VLA observations (iBruni et al.ll2012b . resulting in 
linear sizes of 217, 380, 69, and 17 kpc, respectively. The re- 
maining 84% of the sample turned out to be unresolved, with a 
linear size below 20 kpc in 90-95% of the cases. A previou s 
VLBI study by our group (iMontenegro-Montes et al.ll2008bh . 
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ba sed on a fainter sample of radio-l oud BAL QSOs presented 
in iMontenegro-Montes et alj (l2008al) . shows that on the parsec 
scale a variety of morphologies are found, implying different 
orientations. In that work a total of five sources were observed 
with the VLBA, resulting in two resolved sources (40%). In this 
work we present pc-scale observations of eleven sources: nine 
of them present a resolved structure (~82%) without a predom- 
inant morphology. Both double sources with a radio spectrum 
consistent with being young (and implying jets oriented to large 
angles to the line of sight) and core-jet objects (with their axis 
quite aligned to the line of sight) have been found, reinforcing 
the idea that orientation is not the dominant factor in determin- 
ing the properties of BAL QSOs. In our pc-scale images, three 
sources (0756+37, 1014+05, and 1406+34) can account only 
for part of the VLA flux density measured at the same frequency. 
This missing flux density could be related to some emission with 
angular scales not sampled by the VLBI baselines (i.e. exceed- 
ing several tens of milliarcseconds, and therefore resolved out by 
the observations presented here) but contributing to the arcsec- 
ond scale measurement. Otherwise, such a difference may arise 
from intrinsic variability of the most compact components. All 
of these three sources were unresolved in the VLA maps, with 
linear size below 8 kpc. 

If we also consider the sources 1 159+ 01 and 1624+37, stud- 
ied by iMontenegro-Montes et al.|[2008bl we end up with a to- 
tal of 11 out of 13 (~85%) sources that appear resolved on 
the pc-scale. In 1159+01 a central core was found, with two 
faint symmetric extensions and another two components toward 
SW, while in the case of 1624+37 a core -jet structure was de- 
tected. This is consistent with the variety of morphologies we 
have found in the present work, and it suggests a variety of ori- 
entation with respect to the line of sight, in agreement with the 
results in lBruniet aI1(l2012l) . 

These results, as a whole, agree with the findings available 
in the literature, both in terms of unresolved sources fraction and 
in possible morphologies (and orientations). Figure |4]shows the 
integrated flux density from the FIRST versus redshi ft for the 13 
objects with VLBI observations in the sample from lBruni et al.l 
(2012i) and from the literature. Both redshift and flux density 
ranges are comparable to the literature, so we can safely include 
the results from this paper in the context of previous VLBI ob- 
servations of BAL QSOs. 

After a decade of VLBI observations, BAL QSOs do not 
seem to be predominantly unresolved, or to present any preferred 
morphology. All of them have been detected on the pc-scale. In 
some cases these sources can be compact, with a projected lin- 
ear size below 20 kpc, or even extremely compact, < 1 kpc, but 
there are examples of extended sources (hundreds of kpc), like 
0849+27, presenting BAL features, or cases of BAL QSOs with 
traces of a diffuse emission, probably related to a low-frequency 
component remnant of an earlier active phase. Thus, the inter- 
pretation of BAL QSOs as young sources seems to be overly 
simplistic. On the other hand, the variety of jet axis orientations 
with respect to the line of sight we found does not allow us to 
explain this particular subset of QSOs as likely to be seen from 
a particular line of sight, as suggested by the orientation model 
(Elvis 2000). The outflows at the origin of the BAL features 
should thus be present at various orientations with respect to the 
jet axis, and in various phases of the QSO evolutionary track. 

We determined the magnetic field intensity in the vari- 
ous emitting regions, under the conventional assumptions of 
pure synchrotron emission and minimum total energy (i.e. 
the equipartition field) , and we found values of a few mG. 
iDaUacasa et"an (l2002allbh present VLBA, EVN and Merlin ob- 



servations of a sample of compact steep-spectrum sou rces (CSS) 
drawn from the B3-VLA sample (IVigotti et al.lll989l) and cover- 
ing about the same radio luminosity range of the BAL QSOs 
presented here. They calculated the minimum energy density 
and the magnetic field strength, adopting the same assumptions 
and methodology, for radio components with sizes and flux den- 
sities comparable to the ones presented in this work. They found 
u,„i„ in the range between 10"** and 10"^ erg/cm"*, corresponding 
to B' < 3 mG, comparable with our findings. This suggests that 
the radio-emitting plasma in radio-loud BAL QSOs has physical 
conditions similar to other (non-boosted) sources, while typical 
values for large-scale radio s ources, such as FRI and II, spans 
from tens to hundreds of ixG (jCroston et al.ll2005HStawarz et al.l 



l2005HStawarz & Kataokall2005l) 



5. Conclusions 

We have presented the results of EVN and VLBA observations 
of 1 1 among the brightest radio-loud BAL QSOs in our sample. 
A pc-scale imaging, together with a spectral index analysis of the 
components, was realized. We can summarize our conclusions 
as follows. 

- Nine out of eleven sources (82%) present a resolved struc- 
ture, and various morphologies are visible: double, core-jet 
and symmetric structures have been found, and therefore 
different orientations can be inferred. The percentage of 
resolved sources is comparable to the general results from 
earlier works in the literature. 

- The projected linear sizes of the sources presented here range 
from a few pc (upper limit for the unresolved sources) to 
several tens to several hundred pc. A few sources have addi- 
tional radio emission on much larger scales, as is visible in 
the VLA observations, and therefore a typical radio size of 
BAL QSOs does not exist. These objects explore the same 
parameter space as the common radio sources (ra dio galaxies 
and radio quasars). In the original sample from iBruni et al.l 
(12012), ~80% of sources have linear sizes below 20 kpc, 
still indicating that, on the arcsec-scale, the majority of BAL 
QSOs remain unresolved. 

The missing flux density between VLA and VLBI observa- 
tions can suggest in some cases a diffuse emission that is 
too faint to be detected with these observations. In one case 
(1406+34) this can be due to a low-frequency component, 
probably the remnant of an earlier radio activity of the AGN. 
A complex history can also be supposed for other sources, 
when a re-oriented jet seems to be at the origin of the 
complex morphology found. 

- The variety of linear sizes and spectral characteristics found 
in these VLBI ob servations and from our previous work 
(iBruniet al.ll2012h seems to exclude a simple explanation 
for all BAL QSOs as young compact objects. At the same 
time, given the variety of morphologies that can be found, a 
particular orientation is unlikely to be present: the outflows 
triggering BAL features are then likely to be present in both 
young and older QSOs, at various orientations with respect 
to the jet axis. 

- A calculation of the minimum energy density and magnetic 
field strength shows similar values to the ones present in the 
literature, suggesting that no particular conditions are present 
in the radio-emitting regions of BAL QSOs. 
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